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Evidence for cosmic evolution in the spin of the most 
massive black holes 
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ABSTRACT 

We use results from simulations of the production of magnetohydrodynamic jets 
around black holes to derive the cosmic spin history of the most massive black holes. 
We assume that the efficiency of jet production is a monotonic function of spin a, as 
given by the simulations, and that the accretion flow geometry is similarly thick for 
quasars accreting close to the Eddington ratio and for low-excitation radio galaxies 
accreting at very small Eddington rates. We use the ratio of the comoving densities 
of the jet power and the radiated accretion power associated with supermassive black 
holes with m.^lO 8 M Q to estimate the cosmic history of the characteristic spin a. The 
evolution of this ratio, which increases with decreasing z, is consistent with a picture 
where the z~0 active galactic nuclei have typically higher spins than those at z~2 
(with typical values a~0.35-0.95 and a~0.0-0.25 respectively). We discuss the impli- 
cations in terms of the relative importance of accretion and mergers in the growth of 
supermassive black holes with m.iJfO 8 M . 

Key words: galaxies : active galaxies : jets-galaxies: nuclei - quasars: general -black 
hole physics - cosmology: miscellaneous 



1 INTRODUCTION 

Astrophysical black holes are described by two parameters, 
mass m. and spin a, defining the structure of space-time 
within regions close to the event horizon. Every massive 
galaxy has a supermassive black hole (SMBH) at its centre, 
and the combination of accretion and coalescence of SMBHs 
of similar m ass, major merging, can lead to a wide spread 
in spins (e.g. iBerti fc Voronterill200cf ). A possible constraint 
on the cosmic spin history is the ratio of kinetic to radi- 
ated outputs in active galactic nuclei (AGNs), since these 
are powered by SMBHs. 

During accretion onto black holes, the amount of energy 
available for radiation and for the production of jets is deter- 
mined by the mass and spin of the black hole, as well as the 
rate of accretion of matter and the geometry of the accretion 
flow. More massive black holes can accrete more matter, but 
those that are spinning rapidly can extr act energy more ef- 
ficiently from the infalling mate rial (e.g. iNovikov fc Thorng 
1 19731 : IBlandford fc Znaieklll977h . In addition, geometrically- 
thick accretion flows can power jets more effectively than 
thin ones because they are able to sust ain more powerful 
poloidal mag netic fields (e.g. iMeierlkuOll ). 

In galactic black holes (GBHs), the presence of steady 
jets occurs typically during low accretion rates (where by 
low we mean an Eddington ratio A=m/rhEddSSlO -2 ). Dur- 
ing transitions to moderate accretion rates, AjjlO" 2 the jet 
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is temporarily enhanced in powe r, and can then b ecome 
steady, transient or totally absent (|Fender et alj|2004l b This 
is interpreted as a tr ansition from advecti on-dominated ac- 
cretion flows (ADAFs lNaravan fc Yilll995l ) that are geomet- 
rically thick and can sustain significant jets, to gas-pressure- 
dominated accretion flows which are th in, with aspect ra- 
tios ( the ratio of height to radius) ~0.01 (|Novikov fc Thornd 

Eazjj. 

Quasars and high-excitation radio galaxies are SMBHs 
accreting at a high fraction of their Eddington limiting lu- 
minosity, so that the acc retion flow is geometrically thick 
l|Novikov fc Thornei fl97^ . note GBHs with A~l also pro- 
duce continous jets, e.g. GRS 1915+105, Fender et al. 2004). 
Low-excitation radio galaxies are SMBHs with very low ac- 
cretion rates, ADAFs, that have similarly thick accretion 
flows. The geometry is thus similar at both high and low 
accretion rates. 

SMBHs are known to display a wide range of jet powers 
with no observable differences in the accretion flow, suggest- 
ing that a hidden variable must be controlling the jet power. 
The 'spin paradigm' assumes the black hole spin to be the 
physical pa rameter controlling the kinetic output in the fo rm 
of jets (e.g. IBlandford fc Znaiekll 19771 : ISikora et al.ll2007l) . 

GBHs show a similar wide spread in jet powers, but 
these are not found to correlate with t he published es- 
timates of the spins l|Fender et al.l l2010h . If the spin es- 
timates are correct, this provides very stringent evidence 
against the spin paradigm for GBHs. However, due to 
the uncertainties in the measurements of both the black 
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hole spin and the total jet power, the results found by 
iFender et ail (|201Ch do not provide robust evidence against 
the spin paradigm. The reader is ref erred to Section 7.2 of 
iMartfnez-Sansigre fc Rawlings! |201ll . hereafter MSR11) for 
a detailed discussion. 

In this letter we make the assumption that spin is indeed 
an important factor in powering jets, following the results of 
general-relativistic ma gnetohydrodynamic simulations (e.g. 
lHawlev fc Krolikll2006l ). Hence, assuming the spin paradigm 
to be correct, we present an observational constraint on the 
evolution of the mean spin of the most massive SMBHs as 
inferred from a measured cosmic change in the ratio of the 
power output in jets to the power radiated by active galac- 
tic nuclei (AGNs). We adopt a ACDM cosmology with the 
following parameters: h = _ffo/(100 km s" 1 Mpc -1 ) = 0.7; 
n m = 0.3; tt A = 0.7. 



2 ASSUMPTIONS 

High accretion rates will lead to efficient radiation of the 
thermal energy originating from viscous forces caused by 
magnetohydrodynamic turbulence. A fraction e of the ac- 
creted energy will be radiated away, and the remaining en- 
ergy will be advected into the black hole. The bolometric 
power available for radiation can therefore be described by: 

Lboi = em.c 2 (1) 

where Lboi is the bolometric luminosity due to radiation, 
m. is the rate of accretion of mass onto the SMBH, and the 
term e is the radiative efficiency. 

The shearing of magnetic fields frozen into the accreting 
plasma will allow the extraction of electromagnetic energy 
to power jets. Close to a rotating black hole, the dragging 
of inertial frames will lead to a higher angular velocity (and 
hence larger shear) and additional amplification of the mag- 
netic fields occurs by extraction of energy from the rotating 
black hole. The power available for the production of jets 
can be described as: 

Qjct = r)rh,c 2 (2) 

where Qj e t is the jet p ower and r] is the jet efficiency (e.g. 
lHawlev fc Krolildl200a . this is also the same as tkin in Mer- 
loni & Heinz 2008). Here it is assumed that r\ is a function of 
spin so that rapidly spinning black holes will have a higher 
jet efficiency. 

The AGNs concerned in this letter are powered by 
SMBHs accreting either at very low or at very high fractions 
of th eir Eddington limit i ng luminosity (so either A^10~ 2 or 
A~l IWillott et al.ll200ll : iMcLure et afll2004 ISmolcic et ail 
2009). In these two limits of low- and high-accretion rates the 
geometry of the accretion flows will be similar, with aspect 
ratios <;0.1. Hence, we assume that r\ is a monotonic function 
of a and that variations due to the accretion flow thickness 
are a secondary effect (see also Section 2.2 of MSR11 for 
more details). 

3 THE COMOVING JET POWER AND 
RADIATED ENERGY 

To infer the cosmic spin history of SMBHs, we will use the 
mean comoving kinetic jet power density from AGNs pro- 



ducing jets, PQ jet , and the mean comoving jet power from 
radiating AGNs, pz, bol . These quantities can both be mea- 
sured using: 

PQjet = J Qjot^Q ict (Qjct,z)d<2j ct (3) 
PL bol — J Lbo\4>L hol (Lbo\, z)dLboi, (4) 

where 0g Jet (Qjet, z) = d 2 iV/dQj c tdz and <j>L bol (ibol, z) = 
d 2 A r /dLboidz are the luminosity functions for Qj et and Lboi 
respectively, representing the comoving number density of 
AGNs with a given kinetic jet power or bolometric luminos- 
ity, and their evolution with cosmological redshift z 

The jet power can be estimated from the low-frequency 
radio luminosity de nsity due to synchrotron radiation (e.g. 
IWillott et all 1 1999b . assuming that the jet output results 
in energy stored in radio source lobes together with associ- 
ated and ineluctable work done on the source environment. 
Adopting standard assumptions for the way this energy is 
shared between magnetic fields and particles, the work done 
inflating these lobes as well as typical jet advance speeds 
and a typical density profile for the intergalactic medium, 
PQjct can be estimated from the radio luminosity function 
<j>L v (L v ,z): 

(wltc^)= 3X ^ / G^)V(W>d^ (5) 

where L v is the luminosity density at 151 MHz in 
W Hz -1 sr _1 . The term / represents the combination of 
several uncertainty terms when estimating Qj t from L v (for 
example plasma volume filling factor or the fraction of en- 
ergy in non-radiating particles). An alternative method to 
estimating jet powers is to estimate the energy stored in 
cavities in the X-ray emitting intra-cluster gas , and both 
methods agree best for /~20 (|Cavagnolo et al1l2010l ). We 
therefore initially take /=20 but in Section [4] consider devi- 
ations from this assumption. 

The radio luminosity function includes two components: 
the luminosity func tion for the high rad io luminosity classi- 
cal double sources (jWillott et al.ll200lf ) and the luminosity 
funct ion for the lower radio luminosity AGNs l|Smolcic et al.l 
2009 )Q- For the lower radio luminosity AGNs, the luminosity 
function extends to 2=1.4 only, and we have extrapolated it 
to 2=2.0 using the same functional form for the evolution. 
The luminosity functions of other AGNs suggest a turnover 
around 2=2, so no further extrapolation is attempted. We 
neglect any errors in the quantity inside the integral in Equa- 
tion [5] which is reasonable because it is essentially the lumi- 
nosity density, a quantity dominated by objects near the 
break in the radio luminosity function: t he normalization at 
this break is directly m easured to z ~ 2 (|Willott et al.ll200"ll ; 
IVardoulaki et alj|20ich . 

In AGNs the hard (high-energy) X-ray luminosity orig- 
inates from inverse-Compton scattering of photons from the 
accretion flow by a corona of hot electrons and it is a good 
tracer of the bolometric luminosity. The X-ray luminosity 

1 The latter has been converted from 1.4 GHz to 151 MHz as- 
suming L v ocu~ a with a radio spectral index of a=0.75 
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function 4>l k {Lx, z) 
estimate pL bol ■ 



PL ba 



l|Silverman et al.l 120081 ) can be used 



to 



W Mpc 



= / (1+Fct) Cxix 0L x (ix,2)d£ x . (6) 



The factor Fct accounts for the fraction of luminous 
Compton-thick AGNs missed by the hard X-ray surveys: 
(Fct »0.5, e.g. IMartmez-Sansigre et all 120071 ; ibilli et all 
20071). and we adopt an uncertainty (from Table 8 of 



Fiore et all 120091 ). so that F CT = 0.5 ± 0.1. The bolomet- 
ric correction Cx converts from a monochromatic lumi- 
nosity to a bolometri c luminosity, Lh \ = Cx^x: we use 
the values given by dHopkins et all 120071 ) with, following 
IMartmez-Sansigre fc Tavlorl ( 20091 ). an assumed 8 per cent 
uncertainty. 

The value of PQ jot is always dominated by the output 
from the most massive SMBHs, with m.> 10 8 M Q at all 
redshifts l|McLure et al.l I2004J : ISmolcic et al1l2009l ). To en- 
sure that pL bol is also dominated by the most massive black 
holes, we only integrate above X-ray luminosities that can be 
achieved by SMBHs with m.^10 8 M accreting at >25% of 
the Eddington accretion rate. This is a compromise between 
completeness to SMBHs with m.>10 8 Mq and contamina- 
tion from lower-mass SMBHs. This value securely avoids 
serious contamination from SMBHs with masses <10 8 Mq . 
From the distribution in A found by iHeckman et all <|2004l . 
their Figure 3) we estimate that it leads to pL bol at low-z 
being underestimated by a factor ~3. At high redshift we 
expect it to underestimate pL bol by a negligible amount. 

This evolution of the ratio PQ jot / PL bol with redshift z 
for the most massive SMBHs (with m.il 10 s Mg) is shown 
in Figure 1. There is likely some exaggeration of the rise in 
PQ mt /pL bol towards low z, because of the cut-off assumed in 
X-ray luminosity, but only by a factor ~0.5 (base-10 loga- 
rithm), smaller than the inferred cosmic rise towards low z. 
Note also that the absolute value of the height of the curve 
will scale with / as f 3 ^ 2 . 



4 THE COSMIC SPIN HISTORY 

The ratio of PQ jBt and pL bol , for the most massive black 
holes (<; 10 s Mq) is then given by 



T]Pr, 



(7) 

PL bol e P™.z. bol 

that includes the comoving density of accretion of the AGNs 
contributing to PQ jct , namely Pm, Q . t and of AGNs con- 
tributing to p Lbol , namely Pm. Lbo] ■ 

Separating the contribution to pq . % from ADAFs and 
quasars, pq jot —i]auaf pm ADAF + VQSO prh QSO ■ In the case of 
PL bol , although Pm ADAF ^Pm QSO , £ADAF<£QSO, so that p Lbol 
is totally dominated by quasars and PL bol =£QSOPmQ SO - 

Hence, 



PL b , 



I TM-DAF 
V 6QSO 



+ 



?7QSO 
£QSO 



(8) 



We consider two approaches to using pq . /pL bol to con- 
strain the cosmic spin history. The first one is to attempt to 
use the ratio to find a fiducial spin, under some simplifying 
assumptions. The second is to compare a parametric model 
to the observations, to see how well it reproduces them. 
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Figure 1. The posterior probability distribution function for the 
logarithm (base 10) of PQ jet / PL hal as a function of redshift for the 
most massive black holes (with m m ?J 10 s Mq). The colours rep- 
resent the normalised posterior probability of PQ jct / PL bol , given 
the data. Light yellow represents the highest probabilities, dark 
red the lowest. Contours are also marked as black lines, including 
the line of maximum posterior probability (solid), as well as the 
±1ct (dashed) and ±2<r (dotted) lines. At all redshifts we only in- 
tegrate above X-ray luminosities that can be achieved by SMBHs 
with m.^lO 8 Mq accreting at ^25% of the Eddington accretio n 
rate (the characteristic value found bv lMcLure fc Dunlop 2004). 



In both approaches, the jet efficiency r\ is described 
using the results from a set of three-dimensional general- 
relativistic magne tohydrodynamic simula tions, which can be 
approximated as (|Hawlev fc K rolik 200j]): 



??(|o|) « 0.002(1- |a|) _ 



(9) 



For spin values a— [0.0, 0.5, 0.998] the corresponding jet 
efficiencies are r;=[0.002, 0.004, 1.00]. However, our results 
are very similar if we use other efficiencies from the literature 
(see supplementary figures and MSR11). 

For the ra di ative efficiency, we assume the 
iNovikov fc Thornd l| 19731) model, which for spins a= 
[0.0, 0.5, 0.998] has values e=[0.057, 0.082, 0.321] for 
co-rotating accretioijf]. 



4.1 The fiducial cosmic spin 

At Z/Zl, the density of accretion onto ADAFs is insignificant 
compared to the density of accreti on onto quasars (e.g. se e 
Figure 9 of MSR11 and Figure 4 of iMerloni fc Heinzl [20081 ) . 
PA ADAF «pm QSO so that PQ jct /pL bol =??QSo/eQSO. There- 
fore, at high redshift pg jct /pL bol gives a good estimate of the 



2 In Section [5] we will discuss the relevance of chaotic accretion. 
In this paradigm, an avera ge efficiency of c o- and counter-rotating 
accretion should be used l lKing et al . 2008). The results of MSR11 
used jet efficiencies from the literature for co-rotating accretion, 
so for consistency we limit ourselves to these. The radiative effi- 
ciency is determined for quasars and their spins are close to a=0.0 
(MSR11) so the co- and counter-rotating radiative efficiencies are 
virtually identical. Therefore, although there is a slight concep- 
tual inconsistency between our model radiative efficiency and our 
interpretation, in practice it makes no difference to the results. 
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ratio of efficiencies for quasars and, under the spin paradigm, 
of their typical spin. 

At very low redshifts the comoving density of ac- 
cretion onto ADAFs is comparable to tha t of quasars, 
so /5a adaf « Pm Q so (|Merloni fc Heind 120081 , MSR11), in 
which case PQ jct /pL bol ~ (j?adaf/??qso + 1) ^qso/cqso- It 
is therefore not straightforward to interpret PQ jBt / PL bol 
in terms of typical efficiencies, but some progress can be 
made under the assumption that jjadaf ^QSO, so that 
PQjot/pi-boi^ 2 7?Qso/eQSO- This shows that at z~0 the ra- 
tio PQjct/Piboi wm overestimate the typical efficiencies of 
quasars and ADAFs. 

Figure 2 represents a fiducial cosmic spin history of the 
most massive black holes inferred from PQ jBt / PL bol - It repre- 
sents the characteristic spin if at a given redshift all SMBHs 
have the same spin. It is derived by assuming the first term in 
Equation [8] is negligible, so we expect this figure to provide 
an accurate estimate at moderate and high redshifts, but 
to overpredict cifid(z) at the low redshifts, marked in gray. 
Overestimating t)/e by a factor of >2 at z~0 suggests that 
the typical spin at z~0 is <0.9 rather than ~0.95. Consider- 
ation of / suggests that lower values of / also lead to lower 
estimates of the typical low-z spin: re-running our analysis 
with / = 10 (rather than / = 20) leads to no changes at high 
z but a typical low-z spin of 0.8. We note that high values 
of / are supported by X-ray cavity observations undertaken 
at low radio powers and low redshift so lower values of / 
for the low-z population are hard to envisage. If we assume 
that different effects (e.g. / and underestimating pz, bol at 
low-z) add to an overall overestimate of PQ et / pL bol by a 
factor of ~2 (4), this will result in a fiducial spin at z~0 of 
&fid~0.9 (0.7). In these cases the inferred evolution would 
be less extreme, but still present. 

To test wether this fiducial spin history is con- 
sistent with constraints fr om the radiative efficiency of 
quasars (the ISoItanl 1 19821 . argument), we compute the 
luminosity-weighted spin (see also MSR11): 
j and{z)pL bol dz/ J pL bol dz. The spin history of Fig. f2]yields 
a luminosity-weighted spin is (ci}l=0.055, which corresponds 
to a radiative efficiency of (e)L=0.059. This is in good agree- 
ment with observational const raints, which typically yield 
value s (e)L~0.05-0.10 (see e.g. iMartmez-Sansigre fc Taylor] 
l2009h . 

4.2 Can we explain PQ jot /pL bol ? 

In MSR11, we inferred a moderate evolution amongst the 
SMBHs with m.^10 8 Mg. This evolution was really driven 
by the change in relative space densities of two populations: 
the high Eddington-rate SMBHs, or quasars, which dom- 
inate at high redshifts, were found to have spin distribu- 
tions centred around low spin, a~0. The low Eddington-rate 
SMBHs, or ADAFs, become significant at low redshifts, and 
these were found to have a bimodal spin distribution, with 
a peak centred at <z~0 and another at a~l. These spin dis- 
tributions allowed us to explain the local radio luminosity 
function, as well as explaining its evolution up to 

The results of MSR11 suggest that the evolution of 
the mean spin of SMBHs is driven by a transition from a 
high-accretion rate, low spin population, to a low-accretion 
rate population showing a bimodality in spins (and hence 
a slightly higher mean spin). This evolution was qualita- 




z 



Figure 2. Fiducial spin of the most massive black holes as a 
function of redshift. It shows the posterior probability distribu- 
tion for afl(j(z) given the ratio PQ jBt / PL bal plotted in Figure 1, 
given r?(|a|) and e(|a|). The posterior probability has been nor- 
malised at each redshift. The underlying assumption is that the 
aspect ratio of the accretion flow is similar for all the AGNs con- 
tributing to PQ jet /pL bol ) so that geometry has a negligible effect 
in the evolution of PQ jot / PL bal t which is determined by the evo- 
lution in cosmic spin. The colours and contours are the same as 
those in Figure 1. At z^JO.9, the maximum posterior probability 
line as well as the -1 and -2<r contours go along the a=0 axis. 
The grey shades mark the regions where the assumption that 
PA ADAF < Pm QSO breaks down. In the light gray region, we es- 
timate that Prh ADA F~ 0-lPwiQso while hr the dark gray region, 
P™-adaf~ 0-25prfiQ SO (MSR11). Hence in these regions we over- 
estimate OfldO*)' 



tively similar to that seen in Fig. [2] but it was significantly 
weaker, changing from (a)~0.25 at high redshift to ~0.35 at 
low redshift. 

We can test whether the spin distributions we in- 
ferred in MSR11 can reproduce the observed evolution in 
pQ ct / pL bol - We use these spin distributions to predict how 
the ratio PQ jet / pL bol evolves with redshift using 

( PQjet \ _ f (t/adaf) \ / Paadaf \ + ( (?7QSoA ^ 

XAbol/pred V ( £ QS0) J \ Pm QSO / \(£QSo) J ' 

Note that the angled brackets indicate expectation values: 
(i)adaf) = J ry(a)P A DAF(a)da, (tjqso) = / r?(a)P QS o(a)da 
and (cqso) = / e(a)-PQSo(«)da. The variance of the ADAF 
term is large, due to the bimodal nature, while the variance 
of the QSO term is small. As mentioned earlier, at zjjl 
P m Qso ^P^adaf so that the t]qso term dominates, but at 
low redshift Pm QSO ~Pm ADAF • 

Fig- EJshows the resulting evolution of PQ ct / pL bol , from 
Equation 1101 and using the spin distributions from MSR11. 
At z~0 the most likely predicted value of PQ jct / pL bol is 
about 0.1 dex higher than the observed value, and at zjjl 
it is about 0.3 dex higher. This can be understood since in 
MSR11 (Figs. 4 and 7 of that work) we overpredicted the 
radio luminosity function of high-Eddington rate SMBHs, 
which dominate pQ jct at high redshifts. 

The expectation value is systematically slightly higher 
than the most likely value of PQ ot /pL bol shown in Fig.[T] but 
the curvature is similar and the two figures agree very well 
within their uncertainties. Similar results are found using 
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Figure 3. The evolution of PQ iet / PL boi as predicted by the spin 
distributions from MSR11. In MSR11 we predicted spin distribu- 
tions for high-Eddington rate objects (QSOs) and low-Eddington 
rate ones (ADAFs), as well as the accretion density onto these 
two sub-populations (PriiQ SO and Prh ADAF respectively). Follow- 
ing Equations 151 and 1101 and the spin distributions from MSR11, 
we predict an expectation value that behaves very similarly to 
that of Fig. [T] The contours are wide due to the variance of the 
spin probability distributions for high- and low-Eddington rate 
objects. At low redshifts the variance is larger due to the bimodal- 
ity in the spin distribution of the low-Eddington rate objects. 



other jet efficiencies from the literature (see supplementary 
figures). 

Although the amount of evolution in the mean spin is 
different, Fig. [2] and the mean spin from MSR11, which can 
reproduce Fig. [1] agree qualitatively: the typical spin was 
very low at high redshifts, and has increased at low redshifts. 



5 DISCUSSION 

Using the observed ratio between the comoving densities of 
jet power and radiation from SMBHs, PQ jmt /pL bo} , we have 
inferred that the spin of the SMBHs must be evolving. We 
have done this in two ways. The first method included some 
simplifying assumptions to allow us to "invert" the problem. 
This allowed us to infer a fiducial spin for all SMBHs, and we 
find this fiducial spin evolves strongly. The second method 
consists of comparing some highly parametrised description 
of the spin distributions, inferred from fitting the local radio 
luminosity function in our previous work (MSR11). There 
are no simplifying assumptions, but rather than infer the 
evolution of a fiducial spin value, we simply test whether 
spin distributions previously inferred in MSR11 can explain 
the observed evolution of PQ-. et / PL bol ' we find that they can 
explain this evolution very well. 

The results from both methods are quantitatively differ- 
ent, in that they differ in the spin values inferred at low red- 
shift. However the two methods agree qualitatively: a low- 
spin epoch occurs at z <; 1 and coincides with the period 
during which the black h oles with masses m.^ 10 8 Mq ac- 
creted most of their mass (|Yu fc Tremainell2002l ). This epoch 
is one characterised by high accretion rates onto low-spin ob- 
jects, while the present-day epoch (z~0) is characterised by 



low accretion rates onto a population of SMBHs where a 
significant fraction have high spins. 

Somewhat counter-intuitively, the epoch of lowest-spin 
corresponds to the epoch of highest accretion. This, how- 
ever, is consist ent instead with the paradigm of "chaotic 
accretion" fe.g. lKing et al.ll2008l ): approximately half of the 
accretion occurs with the disc counter-rotating with respect 
to the SMBH, hence decreasing the spin of the hole as op- 
posed to the increase brought by a co-rotating disc. 

Chaotic accr etion is expected to lead to SMBHs with 
low spins, a~0.1 (|King et al.l l2008). and our result is thus in 
agreement with black holes being "spun down" at high red- 
shift, when the typical accretion rate was higher, due to the 
greater available supply of concent rated cold gas in the cen- 
tral r egions of the host galaxy (e.g. lObreschkow fc Rawlingsl 
120091 ). 

When two black holes of similar mass coalesce, the or- 
bital angular momentum contributes significantly to the fi- 
nal angular momentum of the coalesced black hole, lead- 
ing to high values of the final spin, typically a~0.7 
ijRezzolla et al.l 2008). Hence, mergers of SMBHs of similar 
mass (major mergers, with mass ratios 5s4:l) are expected 
to provide the "spin-up" mechanism. 

SMBHs of mass JilO 8 Mq are hosted by galaxies with 
bulge masses ^lO 11 Mq, typically elliptical galaxies. A sig- 
nificant fraction of t hese galaxies have u ndergone major 
mergers since z ~ 1 (jRobaina et al.ll2010l ). meaning their 
central SMBHs will have most likely undergone a major 
merger too. 

Due to the steepness of the mass function above the 
break in the mass function, our result is most likely domi- 
nated by the behaviour of SMBHs hosted by galaxies with 
masses ~1 x 10 11 Mq and reflects their merger history, 
rather than reflecting that of the rarer, more massive galax- 
ies with >3 x 10 Mq, which seem to have undergone an 
earlier evolution |Banerii et al.ll20ld ). 

The mechanism for spinning-up SMBHs, major merg- 
ing, is present during the entire epoch of ^ z ^ 2 since the 
merging of galaxies with similar masses is also a common 
occurrence at redshifts z~2. However, the centrally concen- 
trated, typically molecular, cold gas content of the merging 
galaxies, the fuel for accretion, decreases significantly with 
cosmic time. 

It is therefore the gradual disappearance of the "brak- 
ing" mechanism, chaotic accretion, which determines the in- 
crease in the typical spin of the most massive black holes. 
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